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EXPERIMENTAL AND ANALYTICAL INVESTIGATION OF 


PROPELLER WHIRL FLUTTER OF A POWER PLANT 
ON A FLEXIBLE WING 

By Robert M. Bennett and Samuel R. Bland 
Langley Research Center 


SUMMARY 


Results of an investigation of propeller whirl flutter are presented for a 
model consisting of a single propeller and simulated power plant mounted with 
flexibility in pitch and yaw to a cantilevered; semispan wing. Several config- 
urations differing in the power-plant parameters were investigated analytically 
and experimentally. An analysis is presented employing four uncoupled modes: 
wing bending; wing torsion; power-plant pitch; and power-plant yaw. Experi- 
mental results are compared with this four-mode analysis and with an analysis 
including only power-plant freedoms. 

The experimental and analytical results indicate that the effects of the 
wing on the whirl flutter boundary were large for some eases; depending on the 
system parameters; and were generally stabilizing for the configurations con- 
sidered. The four-mode analysis generally gave better results than the two- 
mode analysis, but only gave fair agreement in some instances. The analyses 
indicated that in cases of large- amplitude wing motion; wing aerodynamics can 
have significant stabilizing effects. 


INTRODUCTION 


A design consideration for propeller-driven aircraft is the prevention of 
a dynamic instability; generally termed propeller whirl flutter or autopreces- 
sion; in which the propeller hub wobbles or executes a whirling motion. The 
flutter for an isolated power plant has been treated analytically in refer- 
ences 1 to 7- Comparisons of measured and calculated flutter boundaries for a 
windmilling propeller (ref. 8) indicated that the whirl-flutter boundary could 
be satisfactorily predicted if accurate aerodynamic parameters were used. Cal- 
culations for thrusting propellers in the cruise condition (refs. 1 to k) have 
indicated little effect of thrust on whirl flutter. However, the question of 
the effects of the wing has been considered only briefly in these investiga- 
tions. Furthermore, comparison of the theoretical results of reference 1 for 
an isolated power plant with the data of reference 9 for a four-engine model 
indicated large differences that might be attributed to the wing and complete 


model. Since further investigation of the effects of the wing was desirable, 
the present investigation was undertaken with the following objectives: 

(1) To determine the whirl-flutter boundaries of a power plant (with wind- 
milling propeller) on a cantilever, semispan wing, which not only has realistic 
parameters but is also more readily amenable to analysis and physical interpre- 
tation than a complicated four-engine model. 

(2) To investigate the system of reference 8 on the semi span wing and eval- 
uate the range of applicability of the theory for an isolated power plant in a 
limited manner. 

( 3 ) To analyze combined wing and power-plant flutter, evaluate the analysis, 
and explore the dynamics of wing — power-plant coupling. 

For this investigation the model of reference 8 was mounted on a canti- 
levered semispan wing having an aspect ratio of 6.97 (based on area of exposed 
wing panel) and a taper ratio of 0.^30- Flutter boundaries were measured over 
a range of power-plant stiffnesses and damping levels and propeller blade 
angles. In addition, three other configurations differing in power-plant param- 
eters were investigated over a range of blade angles for a single value of 
stiffness and of damping. Following the list of symbols in appendix A, the 
equations of motion are developed in appendix B for coupled wing — power-plant 
flutter by using an uncoupled-mode analysis and by omitting aerodynamic coupling 
between the wing and propeller. The results of this analysis are compared with 
the results of the experiment and with the theory for an isolated power plant. 

The effects of gyroscopic coupling on the natural vibration modes are briefly 
considered in appendix C. 


APPARATUS AN D TESTS 


Tunnel 

The investigation was made in the Langley transonic dynamics tunnel which 
is a slotted- throat, variable-pressure, single-return wind tunnel having a test 
section l 6 feet square (with cropped corners). It is capable of operation at 
Mach numbers up to 1.2 and at stagnation pressures from near vacuum to slightly 
above atmospheric. Either air or freon can be used as a test medium. Large 
windows are provided for close, unobstructed viewing of the model. The dynamic 
pressure in the test section can be rapidly reduced when flutter occurs by the 
operation of a bypass valve in a channel which connects the plenum chamber sur- 
rounding the test section to the return passage of the tunnel. 

The present investigation was conducted in air at near atmospheric condi- 
tions and at Mach numbers less than 0 . 3 . 
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Model 


General arr angement and mounting *- The model consisted of a semispan wing 
with a single nacelle-power-plant combination with windmilling propeller. 
Although no attempt was made to scale the model from a specific airplane, the 
model parameters are believed to be realistic since the model was adapted from 
the four-engine model of reference 9- The model was cantilevered from the tun- 
nel wall as shown in an overall view in figure 1. This was accomplished by 
sandwiching the root portion of the wing spar with heavy plates (see fig. 2), 
bolting them together, and clamping the assembly in a remotely operable turn- 
table, which was used to adjust the angle of attack of the model (see section 
entitled "Tests" ). In order to minimize effects of the wall boundary layer and 
to fair around the rigid inboard portion of the wing, a simulated fuselage was 
bolted to the tunnel wall. Sealing around the wing- fuselage juncture was accom- 
plished by means of a flexible foam material. 

Simulated power plant .- Four configurations, differing in the parameters of 
the simulated power plant, were tested and are denoted herein as configurations 
A, B, C, and D. The system consisted of a simulated engine and a propeller 
attached to an aluminum mounting beam through a gimbal with spring-restrained 
pitch and yaw freedoms. A variable- speed motor with an eccentrically mounted 
weight on its shaft formed part of the engine mass and served as a shaker device 
for configurations A and B. The schematic for each configuration is given in 
figure 3- Configuration A was the symmetrical model of reference 8. Configura- 
tion B differed from configuration A only in the location of the lead weights, 
which were moved from the rear of the engine to the mounting beam at the same 
distance from the elastic axis (fig. 3(a)). For configuration C, the simulated 
engine mass was also removed from the gimbal structure and mounted on the 
mounting beam at the same distance from the elastic axis (fig. 3(b)). Configura- 
tion D differed from configuration C in the spring and clip arrangements which 
permitted unequal pitch and 1 yaw stiffnesses. The mass parameters for each con- 
figuration are given in tab^Le I. 

The aluminum beam of this model on which the gimbal was mounted was similar 
to that of reference 9, but was stiffened to raise its frequencies considerably 
above the whirl- flutter frequencies by the addition of L- section doublers 
(fig. 3)- For example, cantilever frequencies of the mounting beam with engine- 
propeller configuration A were 29 cps in pitch and 2k cps in yaw. 

The nacelle was made of balsa construction and was fixed to the wing beam 
with a stiff mounting bracket. Clearance for whirling motions was left between 
the nacelle and spinner. 

The four-blade aluminum propeller was that used for the investigation of 
references 8 and 9- The variation of windmilling advance ratio V/nD with 
reference blade angle is given in figure k (from ref. 8). 

Wing .- The wing, shown in figure 2, consisted of a built-up aluminum spar 
for stiffness and balsa pods for airfoil contour. The gaps between the pods 
were sealed with flexible, foam plastic strips. For analysis purposes, the root 
was considered to be fixed at the bolts through the outboard edge of the sand- 
wiching plates, 1.7 inches inboard of the edge of the first pod. The rounded 
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tip was treated as a streamwise tip ending at the tip of the wing spar (fig. 2). 
The aspect ratio was 6. 97 (based on area of the exposed wing panel) and the taper 
ratio was 0.430. The center line of the wing spar, which was the elastic axis, 
was swept forward 4.9°. The nondimens ional distance from the midchord to the 
elastic axis is shown in figure 5- 

Although the elastic axis was swept forward, the pods were formed as stream- 
wise strips (fig. 2). The mass, static unbalance, and inertia about the elastic 
axis of each strip were measured in the streamwise direction at the spanwise 
position of the center of gravity of the strip, and were considered to be the 
average values over the width of the strip. (The measurements for the spar were 
made on a similar wing spar cut into strips.) The spanwise distributions of 
these parameters are presented in figure 6 in which the nondimens ional parameters 
are based on the streamwise semichord at the spanwise station of the center of 
gravity of the strip. The parameters for the propeller, power plant, nacelle, 
and mounting beam (not included in fig. 6) are presented in table I. 

The spanwise distributions of bending and torsional stiffness are presented 
in figure 7. The stiffness values were obtained by differentiating parabolic 
approximations of the slope of the static deflection curve for a concentrated 
load at the tip. The slopes, which were measured with a mirror system and were 
averaged for several values of loading magnitude, are also presented in figure 7- 


Vibration Modes 

In order to check the mathematical description of the wing by using the 
measured physical properties, a series of vibration measurements were made and 
coit^ared with calculations. For the bare wing spar, which has small static 
unbalance coupling, measured first and second bending frequencies were 17.6 cps 
anrl 59-5 cps, respectively. Corresponding calculated uncoupled frequencies 
(obtained by the method of ref. 10 ) were 17. 0 cps and 66.4 cps. The calculated 
first torsional frequency of the bare wing spar was 267 cps; coupled modes were 
measured at 255 and 270 cps, the latter being more nearly a pure torsion mode. 

With the pods added to the wing spar, but the nacelle omitted, the first 
three measured modes were 10.6 cps (first bending), 39 . 0 cps (second bending), 
and 47.0 cps (first torsion). Corresponding calculated frequencies (considering 
two uncoupled bending modes and one uncoupled torsional mode, all coupled by 
static unbalance in a Rayleigb-Ritz type analysis) were 11.2 cps, 43.1 cps, and 
53.4 cps, respectively. Therefore, the mathematical description of the wing 
using measured physical properties satisfactorily predicts the measured frequen- 
cies of the bare spar or of the wing without the nacelle. 

The measured node lines of the first three modes, which involve primarily 
wing motion, are sketched in figure 8 for the complete model of configuration A 
with f 0 = f^ = 9.3 cps. The calculated uncoupled first bending (10.1 cps) mode 
shape and first torsional (17*1 C P S ) mode shape are presented in figure 9 . When 
these two modes are coupled by static unbalance and the engine pitch mode, the 
calculated coupled frequencies compare with measurements as follows: 
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Predominant 

response 

Calculated coupled 
frequency, cps 

Measured 
frequency, cps 

Engine pitch, 0 

9.2 

9.3 

Engine yaw, i|r 

9-3 

9.3 

First wing tending, h^ 

10.1 

10.4 

Wing torsion, a 

17.4 

14.3 

Second wing bending, hg 

i 


39-1 


Good agreement in frequency is obtained for these modes with the exception of 
the l4.3-cps mode which is primarily wing torsion. The poor agreement for this 
mode may be a result of considering only two wing modes in the calculations. 
Although the next coupled mode involving primarily wing motion has a frequency 
of 39*1 cps, the calculated uncoupled second bending mode frequency is 23.7 cps, 
and its inclusion could possibly affect the results. 


Instrumentation 

An automatic digital readout system was used to record the tunnel static 
and stagnation pressures and stagnation temperature. 

The accelerations of the simulated engine mass were determined by linear 
accelerometers mounted on the nonrotating part of the pivoting mass. The fre- 
quency and a qualitative indication of the amplitudes of the wing bending and 
torsional oscillations were determined by strain gages near the wing root. 

These data were recorded on a direct writing oscillograph recorder for real-time 
monitoring of the test. Propeller rotational speed was measured by a magnetic 
pickup , mounted at the propeller shaft, which drove an electronic counter. 


Tests 

In the tests of configurations A and B, the tunnel was brought up to a low 
speed, the model was initially disturbed sinusoidally by the shaker, and the 
nature of the free vibrations was observed visually and on the oscillograph 
record. Airstream velocity was then increased by a small increment, and the 
model again disturbed. This procedure was continued until a sustained, approxi- 
mately constant amplitude oscillation was produced. The tunnel velocity was 
then rapidly reduced to prevent model damage. The procedure for configura- 
tions C and D was essentially the same; however, since the shakier was ineffec- 
tive for these configurations random tunnel turbulence served as the flutter 
excitation. The absence of the shaker did not appear to impair the selection 
of the flutter condition by the observer. The wing was maintained at nearly 
zero static tip deflection in bending by varying the angle of attack of the wing 
with the mounting turntable. 
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In some cases the stability of the flutter mode was independent of ampli- 
tude as indicated by linear theory. However, as in the tests reported in refer- 
ence 8, cases were encountered in which the amplitude of flutter was limited at 
a given velocity, and increased with increasing velocity, indicating a nonlinear 
effect. The wind-off power-plant damping for these cases increased with ampli- 
tude. The pertinent damping coefficient was assumed herein, as in reference 8, 
to be that derived from the records of the decay of free vibration at the ampli- 
tude of flutter. Also, there were some instances where unstable oscillations 
were produced if sufficiently large amplitudes of motion were excited without 
an increase in velocity. The damping for these cases decreased with increasing 
amplitude} hence, the appropriate damping coefficient was assumed to be that 
derived from the decay records at the amplitude where the motion began to be 
unstable. Damping records were taken before and after each series of runs at a 
given stiffness or damping level. As the power-plant damping varies somewhat 
during a series of runs, the accuracy of 2^ is est i mated to be about ±0.005. 

Power-plant stiffness and damping levels and propeller blade angle were 
varied separately for configuration A. (The damping level was adjusted by 
inserting sponge rubber between the gimbal rings.) Several blade angles, but 
only one level of stiffness and damping, were investigated for configurations B, 
C, and D. A summary of test conditions and data is given in table II. 

All the flutter tests were made in air at near atmospheric conditions and 
at Mach numbers less than 0.5* 


DESCRIPTION OF ANALYSIS 


An analysis of the whirl flutter of a wing— propeller— power-plant combina- 
tion including wing flexibility is developed in detail in appendix B and is 
applied herein to the model of this report. The analysis is based on a 
Lagrangian formulation of the equations of motion employing deflections in four 
uncoupled vibration modes as generalized coordinates. The system treated analyt- 
ically is depicted in figure 10. It consists of a single, rigid engine-propeller 
system mounted with spring-restrained gimbals to a rigid mounting beam which is 
mounted on the flexible wing. Two vibration modes of the wing, bending and tor- 
sion, are considered. These are described as bending and twisting about a 
straight, unswept elastic axis. Both propeller and wing aerodynamic forces are 
included} however, aerodynamic interference effects between the wing and pro- 
peller, and nacelle aerodynamic forces are neglected. The aerodynamic forces 
acting on the propeller are expressed in terms of stability derivatives as in 
references 1, 2, and 8. The modified- strip analysis method of reference 11 is 
employed to represent the wing aerodynamic forces. 

From this analysis the speed for neutral stability (boundary between stable 
and -unstable motion) is determined. This flutter speed, expressed in nondimen- 
sional form as V/Rdq, is determined as a function of the damping of the engine 
pitch mode, 2 £9, by specifying the following damping parameters: g a , g^, and 

2^/2 £q. The method of solution and input quantities required are discussed in 
appendix B. A discussion of the determination of the input quantities for the 
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model of this report is also given in appendix B. The comparison of the results 
of the application of this analysis with experiment for the configuration of this 
report is discussed in the following section. 


RESUITS AND DISCUSSION 


Remarks on Wing Flutter 

A conventional wing flutter .analysis as obtained from the flutter determi- 
nant (eq. (b 43)) by assuming complete rigidity of the power-plant mounting in 
pitch and yaw and by using the first torsional and bending modes of the wing 
indicated no flutter solution with or without propeller aerodynamics. It should 
be pointed out that studies of references 12 and 13 indicate that, in general, 
more than two uncoupled modes would be required to predict flutter accurately 
when large concentrated weights are attached to the structure. Furthermore, 
flexibility in the mounting of a concentrated mass can have an important influ- 
ence on flutter (ref. lb). However, since the concentrated nacelle mass is 
located well ahead of the elastic axis - a condition which has been shown to 
have a strong stabilizing influence on wing flutter (see ref. 13 ) - it is felt 
that the present predictions based on a two-mode analysis are at least indicative 
that classical wing flutter occurs at velocities well above those considered 
herein. Hence, the primary emphasis in this investigation will be on the effects 
of the wing on propeller whirl flutter. 


Whirl Flutter 


The experimental flutter data for configuration A are presented in figure 11 
in the form of viscous damping of the power-plant modes required for stability 
as a function of reduced velocity V/Ro 5 in terms of the average wind-off fre- 

p. 03q + 0Ui f 

quency, ao = — the stable region is above the boundary. Such a boundary 


is unique for a given blade angle for flutter of a power plant flexibly mounted 
on a rigid structure, but becomes a function of the mass, natural frequencies, 
and damping of the wing for the power plant mounted on a flexible wing. Thus, 
the power-plant stiffness level is indicated (fig. ll) by the ratio of uncoupled 
frequencies oS/cr^. As density varied slightly from one test point to another and 


the damping in pitch and damping in yaw were not widely different, the damping 
*2 2£q 4* 2(^k p 0 

parameter 2£ = — has been used to adjust the experimental data to a 


reference density value of 0.0022 slug/cu ft. This relation is an approximate 
one for the power plant mounted on a rigid structure (ref. 8). The power-plant 
damping is treated as viscous damping, as better agreement was obtained in refer- 
ence 8 with viscous damping than with structural damping. However, the damping 
of the wing modes is considered to be structural damping herein. 


The data for configuration A are compared with the theory including only 
power-plant pitch and yaw freedoms (from ref. 8) and with the theory including 
four modes: power-plant pitch and yaw, wing bending, and wing torsion for a 
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power-plant restraint stiffness resulting in a frequency ratio of I.13. 

A comparison of the data with the theory for two modes shows good agreement for 
some cases, hut only fair for others. It should he noted that the results of 
reference 8 for an isolated power plant with two degrees of freedom indicated 
very good agreement between theory and experiment. The theory for four modes 
(for &/o>k = I.13) shows a stabilizing effect in comparison with the theory for 
two modes at the lower damping levels and is in very good agreement with the 
data (for gS/cd^ = I.13) except at a blade angle of 58°- Appendix C indicates 
that the regions of large coupling between the whirling modes and the wing modes, 
as indicated by frequency shifts, are near the coincidences of the frequencies 
of the whirling modes with those of the coupled modes involving primarily wing 
motion. The whirl flutter frequencies for configuration A (table II) were lower 
than those of the first coupled mode; however, for values of to above the first 
bending frequency the whirl flutter frequencies were as much as 0-9 of the fre- 
quency of the first coupled mode. No consistent effect of the higher whirling 
frequencies is shown in the data (fig. 11 ) . 

The flutter boundaries calculated for configurations B, C, and D including 
only power-plant pitch and yaw modes are compared with the boundaries calculated 
including wing bending and torsional modes for three blade angles in figure 12. 
For configuration D, the presentation is based on v/Rdq as cdq and city are 
widely different (table II) . A small stabilizing effect of the wing is pre- 
dicted for configurations B and D, but a large stabilizing effect is predicted 
for conf iguration C at the higher blade angles. Also, for configuration C where 
the wing is signif icantly stabilizing, the flutter boundary is less sensitive 
to damping (fig. 12(b)). 

The theoretical flutter speed boundaries employing either two or four modes 
are compared with experimental data for a single value of power-plant restraint 
stiffness and damping for each configuration in figure 13. Corresponding values 
of flutter- frequency ratio are also presented in figure 13. The theoretical 
values for configuration A, taken from figure 11 , show little effect of the wing 
on the flutter speeds (fig. 13(a)); however, some effect is indicated by the 
experimental flutter speed at a blade angle of 58° • Calculated flutter frequen- 
cies for two- and four-mode analyses differ little and are in good agreement 
with experiment (fig. 13(a)). For configuration B, the theoretical flutter 
speed boundary employing four modes is higher than that obtained using two modes, 
but, as for configuration A, is in good agreement with experiment only at the 
lower blade angles (fig. 13(b)). Calculated flutter frequencies show little 
effect of the wing, but are higher than experiment by about 10 to 15 percent 
(fig. 13(b)). For configuration C, the theoretical flutter boundary employing 
four modes is considerably higher than that calculated with two modes and is in 
good agreement with experiment over the range of blade angles investigated 
experimentally (fig. 13(c)). Flutter frequencies are also in better agreement 
from the four -mode analysis (fig. 13(c)), but are high by about 10 to 25 percent. 
Only a small stabilizing effect of the wing is indicated by the four-mode anal- 
ysis for configuration D (fig. 13 (d)). However, both theoretical boundaries are 
higher than experiment. Calculated flutter frequencies are in good agreement 
for the two-mode calculations, but are low when four modes are employed 
(fig. 13(d)). 
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Thus, for configuration A the wing had a small effectj however, for con- 
figurations B and C the wing had a significant stabilizing effect. The addition 
of the two wing modes to the theory of reference 1 generally improved the agree- 
ment of the theory with experiment, particularly for configuration C, but did 
not fully predict the stabilizing effects or trends. The results of reference 12 
showed that the flutter mode shape of a wing with a large concentrated mass could 
be quite different from the first bending or torsional mode. In particular, the 
relative amplitude at flutter at the position of the concentrated mass was higher 
than that for the first bending mode and lower than that for the first torsional 
mode. Such a result could lead to the observed trends herein, and suggests that 
inclusion of additional wing modes may be required for better correlation. 

In reference 1 a comparison of the theory for the isolated power plant with 
data of reference 9 indicated the possibility of large stabilizing effects of 
the wing and complete model. Also, unpublished analytical studies by R, J. Zwaan 
and H. Bergh of the National Aero-Astronautical Research Institute, N.L.R., 
Amsterdam, .show that generally the wing has a stabilizing effect on whirl flut- 
ter, but under certain conditions wing flexibility can give a reduction in the 
flutter speed. 

Although the theory including four modes was generally conservative for 
configurations A, B, and C, for configuration D, the four-mode theory predicted 
higher speeds than the two-mode theory and was higher than experiment for uncon- 
servative). For configuration D, the ratio of yaw stiffness to pitch stiffness 
was approximately 9 : 1* However, if the flexibility of the nacelle mounting 
beam, which was not considered in the analysis, were taken into account this 
ratio would be reduced somewhat; this may possibly account for the measured 
flutter speeds being somewhat lower than the predicted flutter speeds. 


The Effect of Wing Aerodynamics 

The effect of wing aerodynamics has been considered by omitting the wing 
aerodynamic terms, equation (B39)> In a four -mode analysis for configuration C. 
The resulting flutter characteristics are compared with the two-degree-of- freedom 
calculations and four- degree -of- freedom calculations including wing aerodynamics 
in figure l4 for three propeller blade angles. The asymptotes (taken from 
fig. 15 ) which the coupled frequencies approach for large values of H/cE (with 
no wing or propeller aerodynamics) are included for comparison. These asymptotes 
represent boundaries which the coupled frequencies cannot cross (see appendix C). 
For the 25° blade angle, which is for large values of and beyond the 

region of high coupling, there is little difference among the calculated trends. 
However, for the other blade angles, two flutter solutions were obtained when 
wing aerodynamics were omitted. These two solutions were obtained in the region 
of high coupling, as indicated by the approach of two coupled frequencies to one 
another. However, when wing aerodynamics are included (fig. l4), only one of 
the flutter solutions is obtained. Thus, it is indicated that the inclusion of 
wing aerodynamics is quite important for cases that involve large coupling of 
wing and propeller motion; the omission of wing aerodynamics could lead to an 
unnecessarily conservative flutter speed. 
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CONCLUSIONS 


Propeller whirl flutter boundaries were measured for a single power plant 
mounted on a cantilevered semispan wing for several configurations differing in 
power-plant parameters. An uncoupled-mode analysis of the system was made 
employing the following four degrees of freedom: wing bending, wing torsion, 

power-plant pitch, and power-plant yaw. Comparisons of experimental data, 
results of the four-mode analysis, and the results of an analysis including 
only power-plant freedoms indicate the following conclusions: 

1. The effect of the wing on the whirl flutter boundary can be large 
depending on the system parameters. These effects were generally stabilizing 
for the configurations considered. 

2. The four-mode analysis generally gave better results than the two-mode 
analysis, but only gave fair agreement in some cases. 

3- In cases of large coupling between the wing and the propeller whirl 
modes, the inclusion of wing aerodynamics in the analysis can have large sta- 
bilizing effects. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 28, 1964. 



APPENDIX A 


SYMBOLS 


flutter determinant element (see eq. (B^2) ) 

nondimensional distance from wing midchord to elastic axis, 
measured perpendicular to elastic axis, positive rearward, 
fraction of semichord b 

nondimensional distance from wing midchord to local aerodynamic 
center (for steady flow) , measured perpendicular to elastic axis 
positive rearward, fraction of semichord b 

distance from wing midchord to elastic axis at propeller station, 
positive rearward, ft 

ratio of wing local semichord to reference semichord, measured 
perpendicular to elastic axis, b/b r 

wing semichord, measured perpendicular to elastic axis, ft 

wing reference semichord, O. 63 W- ft 

My p 

propeller pitching-moment coefficient, - 

pV 2 S T R 

M Z p 

propeller yawing-moment coefficient, — ^ — 

pV^S T R 


Fy p 

propeller side-force coefficient, 

ipV 2 S ' 

2 ' 

F Z p 

propeller vertical-force coefficient, - — 

ipV 2 S’ 

2 

wing local lift-curve slope for a streamwise section in steady 
flow, per radian 

dissipation function (see eq. (B2l)); also propeller diameter, 2R 
ft 

wing bending stiffness, lb-ft 2 



e 


a 


e e 

e , 
t 


F 

F Y, P 

F Z, P 
f 

f i 

G 

GJ 

Si 

h 

la p 

*a,W 

^P 

^WP 

I 0,P 


distance from wing elastic axis to plane of propeller-blade 

quarter- chord points at three-quarter radius, positive rearward, 
e 0 + z 0 + (atOp, ft 

distance from gimbal pitch axis to plane of propeller-blade 

quarter-chord points at three-quarter radius, positive rearward, 
ft 

distance from gimbal yaw axis to plane of propeller-blade quarter- 
chord points at three-quarter radius, positive rearward, ft 

real part of complex circulation function based on k r 
propeller aerodynamic force along Y-axis, lb 
propeller aerodynamic force along Z-axis, lb 

vibration frequency, cps 
vibration frequency in ith mode, cps 

imaginary part of complex circulation function based on k r 
wing torsional stiffness, lb-ft 2 

structural damping coefficient for ith vibration mode 

local vertical bending displacement of wing elastic axis, positive 
down, ft 

vertical bending displacement h of wing elastic axis at the 
propeller station, ft 

mass moment of inertia of propeller about axis of rotation, 
slug- ft 2 

generalized mass moment of inertia of wing about elastic axis 
(see eq. (b 6) ) , slug-ft 2 

generalized mass moment of inertia of engine-propeller about 
elastic axis (see eq. (B12) ) , slug-ft 2 

generalized total mass moment of inertia about elastic axis, 

I a,W + *a,P> slug-ft 2 

generalized mass moment of inertia of engine-propeller about gimbal 
pitch axis (see eq. (B12)), slug-ft^ 
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1^ p generalized mass moment of inertia of .engine -propeller about 

gimbal yaw axis (see eq. (B12)), slug-ft^ 

I 0a P generalized cross mass moment-of-inertia term (see eq. (B12) ) , 

9 slug-ft^ 

I 0 *,I. generalized mass moment of inertia of propeller about rotation 

axis (see eq. (Bl6)), slug-ft 2 

i = n 

ig mass moment of inertia of engine -propeller about gimbal pitch axis 

J (see eq.. (B13))> slug-ft^ 

i.. mass moment of inertia of engine -propeller about gimbal yaw axis 

(see eq. (B13) ) , slug-ft^ 

i a w (n) wing mass moment of inertia about elastic axis per unit length 
(see eq. (B7) ) , slug-ft 

J windmilling propeller advance ratio, V/nD 

dimensional stiffness in ith vibration mode 

k r reduced frequency based on wing reference semichord, b^/v 


M p 

% 


distance from wing midchord to gimbal pitch axis, positive 
rearward, ft 

distance from wing midchord to gimbal yaw axis, positive rearward, 
ft 

engine -propeller generalized mass (see eq. (B12) ) , slugs 
wing generalized mass (see eq. (b 6)), slugs 


M WP 

M ij 

My,p 

M Z,P 

Ma(n) 


total generalized mass + Mp, slugs 

nondime nsional generalized mass (see eq. (B27)) 

propeller aerodynamic moment about Y-axis, ft-lb 

propeller aerodynamic moment about Z-axis, ft-lb 

oscillatory wing aerodynamic moment about elastic axis per unit 
length, positive leading edge up, lb 

engine -propeller mass (see eq. (B13)), slugs 
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m p (x) 

ntyU) 

m w (x,y) 

n 

P(ti) 

% 

Q i,P 

Qi,W 

Qij,P 

Q ij,W 

a 

*1 

R 


S’ 

s a,P 

s a,W 

s a,WP 


engine -propeller mass per unit length, slugs/ft 
wing mass per unit length (see eq. (B7) ) , slugs/ft 
wing mass per unit area, slugs/sq ft 
propeller rotational speed, rps 

oscillatory wing aerodynamic lift per unit length, positive down, 
lh/ft 

generalized force in ith vibration mode 

generalized propeller aerodynamic force in ith vibration mode 

generalized wing aerodynamic force in ith vibration mode 

nondimensional generalized propeller aerodynamic force (see 

eq. (B36)) 

nondimensional generalized wing aerodynamic force (see eq. (B39)) 
propeller pitch rate, 0R/V 

ith generalized coordinate (see eqs. (b 4) and (BIO)) 
complex amplitude of q^ (see eq. (B24) ) 

propeller radius, 0.8438 ft 
propeller yaw rate, ijrR/V 

nondimensional radius of gyration about wing elastic axis. 



propeller disk area, itR^ sq ft 

engine -propeller generalized mass unbalance about elastic axis 
(see eq. (B12)), slug-ft 

wing generalized mass unbalance about elastic axis (see eq. (b 6)), 
slug-ft 

total generalized mass unbalance about elastic axis, S a ^ w + S a ^p, 
slug-ft 
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S 0 J ,P 


s a,w(^ 


T 

T P 

T W 

T n 


t 

u 


V 

X,Y,Z 

x,y,z 


y P 

Zp 

Z W 

a 

Op 

p o . 75 R 

S 


engine -propeller generalized mass unbalance about pitch axis 
(see eq. (B12)), slug-ft 

wing semi span , ft 

engine-propeller mass unbalance about pitch axis (see eq. (B13)), 
slug-ft 

wing mass unbalance about elastic axis per unit length (see 
eq. (B7)), slugs 

kinetic energy, ft-lb 

kinetic energy of engine -propeller system (see eq. (b 8 )), ft-lb 
kinetic energy of wing (see eq. (B 5 ))j ft-lb 

kinetic energy of propeller rotation, including gyroscopic pre- 
cession energy (see eq. (Bl^)), ft-lb 

time, sec 

potential energy, ft-lb 

free-stream velocity, ft /sec 

coordinate axis system (see figs. 10 and 1 6 ) 

coordinate distances in X,Y,Z coordinate system, ft 

nondimensional distance from wing elastic axis to local center of 
gravity, measured perpendicular to elastic axis, positive rear- 
ward, fraction of semichord b 

lateral displacement of mass element of engine -propeller system, 
positive right, ft 

vertical displacement of mass element of engine-propeller system, 
positive down, ft 

vertical displacement of mass element of wing, positive down, ft 

local wing angle of attack, positive leading edge up, radians 
wing angle of attack a, at propeller station, radians 
propeller blade angle at three-quarter propeller radius, deg 

virtual displacement in ith generalized coordinate 
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virtual work, ft-lb 


viscous-damping coefficient, ratio of viscous damping to critical 
damping for ith vibration mode 


average viscous-damping coefficient. 


^9 + Pq 

2 P 


nondimen sional spanwise variable, fraction of semispan s 

pitch angle of propeller- shaft axis about gimbal axis, measured 
from wing chord line at propeller station, positive propeller up, 
radians 

effective engine-propeller pitch angle, measured from relative wind 
direction (see eq. (B33))j radians 

jtpR 5 b~ 

density- inertia ratio, — 

I 0,P 


sweep angle of wing elastic axis, deg 
frequency ratio, (JDq/cjd 
air density, slugs /cu ft 

$ 

nominal air density, 0.0022 slug/cu ft 
vibration mode shape in ith mode 


vibration mode shape in ith wing mode at propeller station 


yaw angle of propeller -shaft axis, radians 

effective engine -propeller yaw angle, measured from relative wind 
direction (see eq. (B33) ) ; radians 

propeller rotational speed, 2jm, radians/sec 

vibration frequency, radians/sec 

vibration frequency of ith vibration mode, radians/sec 

cDq + an. 

average wind-off power-plant vibration frequency, - — 1, 

radians/sec 



Subscripts: 

f condition at flutter 

1, j modal designation; i,j = 0, i| f, h, or a 

Partial derivatives are denoted by double subscripts; for example 
oC m dC m 

c mQ = = - — > and so forth. 

Dots over symbols indicate derivatives witb respect to time. 


APPENDIX B 


ANALYSIS 


The analytical model considered here is pictured in figures 10 and l6. 

The system treated has four degrees of freedom: pitch and yaw of a single 

rigid engine-propeller system and wing bending and torsion. A modal-type anal- 
ysis is employed to obtain stability boundaries for the system. The two elas- 
tic modes of the wing are described as bending and twisting about a straight; 
unswept elastic axis. Propeller and wing aerodynamic forces are included; how- 
ever; the aerodynamic interference effects between the wing and propeller and 
the nacelle aerodynamic forces are neglected. 


Equations of Motion 

A Lagrangian formulation of the equations of motion is employed in which 
deflections in the four uncoupled vibration modes serve as generalized 
coordinates. 

Kinetic energy .- The total' kinetic energy of the system is given by 

T = T w + T p + T n (Bl) 

where T w represents the contribution of the wing and rigid power-plant sup- 
porting structure; Tp represents the contribution of the flexibly mounted; 
nonrotating propeller-engine system; T^ is the energy of the rotating propel- 
ler; including the kinetic energy of gyroscopic precession. 

The kinetic energy of the wing is 

T w = I [[ a* &y (B2) 

^ ^ Wing 

where; for bending of and twisting about an elastic axis; 

z(x;y;t) = h + (x - ab)a (B3) 


The bending h and torsion a deflections are given in terms of mode shapes 
and generalized coordinates by 



b(y,t) = <P h (y)q h (t) 


a 


(y.t) = ap a (y) q^Ct) 


(b4) 


Substituting the expressions in equations (B3) and (b 4) into equation (B2) 
leads to the result 


T W - I + § I a ^ 


(B5) 


where 


% = s J dT l 

S a,W = s f Q dT l 

X a,W = S Vw (, l )(f a dT1 






are generalized mass parameters and 


“uU) = / nyCx^yidx 

J Chord 




3 a,w(n) = / rn^(x^y) (x - ab)dx 

J Chord 


i a ,tf(l) = / ny(x,y)(x - ah) 2 dx 


Chord 

are wing- section mass parameters. 

The kinetic energy of the nonrotating propeller is 


J 


Tp = 2 S m p(x)(z| + yp)dx 


(B6) 


(bt) 


(B8) 


where the mass of the system is assumed to be concentrated along the propeller 
axis and the integration is taken over the pivoting system, and where 
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Zp(x,t) = hp + (x - ab)ap + ^x - Zgj( 
y p (x,t) = -(x - Z^t 


The pitch and yaw deflections are given in terms of generalized coordinates by 


e(t) = 9 0 <l 0 (t) 


+(t) = qyi^(t) 


Substituting these results into equation (B8) leads to the result 


1 w *2 


D P = 2 M P q h + S a,PV*tx + 2 + S e,P q 9 q h + Vp^Si 

. 1 T -2.it ’2 

2 I e,p q e 2 ^P^t 


where 


M P = mpcp^p 


D,P = j^e + £e - (ab)pjm p j cp^p^p 

l,p = fe + 2 [ z e - ( ab )p] s e + [ z e - ( ab )p] m vj 


s e,p " ^e^h.p 


|i e + [»o - ( at )p] E eJ 


E e V«,p 


T . m 2 

^P “ 


are generalized mass parameters and 


it in « in 1 1 i ii ■ ill ■■ 



(B13) 


nip 

s e 

i e 

H 



mp(x)dx 




f 

I 

f 


mp(x) (x 

m p( x ) ( x 

mp(x) (x 


Z Q )dx 
hf dx 

|r ) 2 dx 




are engine-propeller mass parameters. The integrations are carried out over 
the pivoting engine-propeller system. 

The kinetic energy of the rotating propeller is given hy an obvious exten- 
sion of the result in reference 1 as 





(B1*0 


In terms of the generalized coordinates in equations (B4) and (BIO) equa- 
tion (Bl4) becomes 


= \ + + ( B 1 5 ) 

where is the mass moment of inertia of the propeller about its axis of 

rotation and 


I i|ra = ^V^P 


(B16) 


are generalized mass parameters of the engine -propeller system. 
The total kinetic energy is seen to be 

T = 2 + 1 0^‘^6 <1 t + S 0,P^0 ( ih + I 0a,P , i0ix + \ I i|f,P^| + 


+ | + s a,wp4 h qa + | i a ,wp^ + | W 2 


(B17) 
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where 


%P - % + M p 


"\ 


3 ct,WP = S a,W + S a,P > 


^(XjWP + lx,,P J 


(B18) 


Potential energy .- The potential energy may be expressed in the form 


° - \ + 1 v4 + 1 + 1 


(B19) 


where K is generalized stiffness. In terms of uncoupled frequencies and gen- 
eralized masses the potential energy is expressed as 


IT = — 


o o 


Dissipation function .- The structural damping concept, in which a dissipa- 
tive force proportional to displacement is employed, will be used in the deriva- 
tion herein. For such damping, the dissipation function is 


2 

“0 -2 


2 

.2 


,2 


CD 


D - | I 0,PS0 W q 0 + 2 I DP S t of ^ + | %P s h a? ^h + | I a,WP g a 


(B21) 


where cd is the frequency of vibration of the system and g is the structural 
damping coefficient. 

Lagrange 1 s equations . - The equations of motion of the system are obtained 
'from Lagrange T s equation 

d /dr \ _ St_ + du_ + 5 d_ _ 
dt^4ij bq_ ± bq_ ± 

where i = 0, i|f, h, or a. Substituting the expressions from equations (B17) , 
(B20), and (B2l) into equation (B22) and performing the indicated operations 
yields the equations of motion 


(B22) 


22 



. p \ 

I 0,py q -e + s 0 oT ^9 + ^ 0 %) + I 0i)r n<1 t + S 0 J P q h + = % 

- I 0 1 )|flq e + ^ - I \(fa fi( ^a = 

s 0,p'4e + M wp(^h + Sh ^ 4 h + + S a,Wp'4x = %. 

I 0a,p4 e + I^4 t + s a,wp'4h + I a,Wp(4a + 8a ^ + <&a) = 






(B23) 


J 


Assume simple harmonic motion with 


ie - So' 1 ”* 

q + = V 11 * 

% - %<=“* 

qa ■ v** 




(B 2 k) 


where q^ , etc., are complex amplitudes of motion. Making these substitutions 
in equation (B23) and dividing the 0 , \|r, and a equations by 

and the h equation by - | I 0 p o) 2 e ia:rt (R being the propeller radius) yields 
the nondimensional equations of motion 


[i - a 2 (i + ig e |5e - iM e+ + M ei ^ + e- 10 * 




iM, 


0 f 


n 4 q + M., 


JTCj. 


l e - «T|n|r 


*T 

v 2 


■ 0 ,P 

1 - X2 (^)V + ig t)J + “* - - 1^? e ' 1Mt 


*eb«6 + “hi 1 - (l + lg h) 


^ + «ha 5 a = ' q- 1 "* 


Mq a 4 0 - i%a jfr 4* + M ha ^ 


I 0 ,I ^ u2 

1 - A 2 (^\ (l + ig a ) 
\cjd q V a / 


a = e " icut l 

^ 

e,P j 

(B25) 


23 



where 


A 


CDq 

CD 


J = 


JtV 



~\ 


(B2 6 ) 


are nondimens ional frequency parameters, h r is the wing reference semichord, 
and 



M 


9h 


s e,p R 

I e,p 


M 


0a 


I 9a,P 

I e,p 




I e,p 


M *a 



M hh 


%p r2 

I 0,P 




s a,WP R 

I 0,P 


Marc 


fa^WP 

I e,p 


-\ 




j 


(B27) 


are nondimen sional generalized masses. 
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Generalized Aerodynamic Forces 

The aerodynamic forces acting on the propeller are expressed in terms of 
propeller stability derivatives as was done in references 1, 2 , and 8. The 
modified- strip analysis method of reference 11 is employed to represent the 
wing aerodynamic forces. The directions of the aerodynamic forces and moments 
are shown in figure l6. 

Virtual work .- The expression for the work done on the system by virtual 
displacements of the coordinates is 


5W = Q e 5q^ + + Q h Sq h + 


(B28) 


where the generalized forces have the following representations in terms of pro- 
peller and wing aerodynamics 


% “ %,T 




"\ 


Q h " Q h,P + Q h,W 


“ Sx,,P + ( ^<x,W 


(B29) 


J 


Propeller aerodynamic forces.- The aerodynamic forces acting on the propel- 


ler are (see fig 


sroayn 

. 1 6 ) 


Q e,p = ( m y,p + e e F z,p 

Q t,P = (^P " e i|r F Y,p)^ 

Q h,P = ^P^P 

Qa,P = ( m Y,P + e aFz,,?)Va,Y 




where 


J 


(B30) 


e a = e 0 + l Q + ( ab )p 


(B31) 
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As in references 1 and 8, My p and Mg p are moments about the Y~ and Z-axes 
in the propeller plane and Fy p and Fg p are forces acting in the Y and 

Z directions on the propeller axis in the propeller plane. These forces and 
moments may be expressed in terms of propeller stability derivatives as 


- 1 ^v2 


■Z,P - 2 - “'(v *\V e * + % I *) 

F Y,P - | 0 v2s ’(<V + % | 5 + %+ + °Y r I t) 


“y,P = PV 2 3'p(c Be 5 + cv lq £ 9 + C„ 1# t 


s-f ♦) 


~\ 


> (B32) 


= pV^S’R^ nQ 0 + £ 0 + Cn^t + C 


R « 




where S' = rfRr , and 


0 = 0 + c&p + 9 + — <ip + i hp 

^ V V ^ V p 


— • 

i(r = ijf + — - ijf 

V 


J 


J 


(B53) 


are pitch and yaw angles with respect to the relative wind direction. Notice 
that in equation (B32) all possible force and moment derivatives that depend on 
angle and angular rate are included; angular acceleration derivatives (which 
account for virtual mass effects) have been neglected, as in reference 1. The 
substitution of equation (B33) into equation (B32), the use of generalized 
coordinates (eqs. (B4) and (BIO)), and the assumption of harmonic motion 
(eq. (B24-) ) lead to the final expressions for the nondimensional propeller 
aerodynamic forces 


^8,>P e -iait 

I e,i^ 2 

e " ia3t 

RQ h,P e _icot 

I 0,I^ 2 

^a,P e -iait 


(%0,P^0 + Q 0\(f,PSf + Q 0h,P ^ + Q 0a,P^aj 

^( Q i|r0,P^0 + Q i|nJf,P^\(r + Q tb,P ~t + Q i(fa,P^a) 

“|( Q h0,P^0 + Q ht,P^|r + Q hh,P + Q ha,p5a) 

~(^a0,P^0 + ^ai|f,P^t + ^ah,P + Sia,P^a j 






J 


(B3^) 


2 6 



where 


Kp 


3 2 
jrpR b r 

I e J ,p 


is an inertia ratio parameter and 


(B35) 
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Wing aerodynamic forces,- The aerodynamic forces acting on the wing 
(introduced in eg. (B 29 )) are given in general as 


Qh,w = s r p ( T i) < Ph. ^ 
d 0 


Q 


a,W 


s 



~\ 




J 


(B37) 


where P(t]) is the wing-section lift force and M^tj) is the section pitching 
moment about the elastic axis. Any suitable unsteady aerodynamic theory may be 
used to evaluate these expressions. In this report, the modified- strip analysis 
method of reference 11 is employed as follows: 


RQ h,W 


Qg„W 

I 0,^ u2 


-iO)t _ L , Q r ' 

2 r*hh,W R ^ha,W^a 




-ia)t 


% 


r \ 


6 = ^l^ah.W f + Qcxa.vAxj 


(B38) 


where, for an unswept elastic axis, 






%h, W = 


G 

jtk r 


+ i 


Jtk r 


sh r 2 

P -1 * 





(B39) 


The circulation functions for unsteady flow, F and G, depend, in general, 
on kj. and Mach number. The circled numbers are the integrals given in 
appendix A of reference 11. They are, in the present notation. 
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1 


(D 


CD 


© 

@ 

© 

© 

© 

© 




b 2 ^ d n 


c z a Bc Ph dT i 


C I B'Ph'Pa dl l 

CL 


2 °*a 


c * a B liT + &C ' *Wa ^ 


B5aCp h. Cp a dT1 


B2 cp h cp a dT! 







J 


(B^O) 


Flutter Determinant 

Substituting equations (B34) and (B 38 ) into equation (B25) yields the 
final form for the equations of motion. In matrix notation 
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where 


CD 

CD 

CD 

< 

A 

CD 

< 

> 

CD 

P 

1 


rv 

A ie 

A \|n|r 

A i|fh 

A \(ra 



CD 

4 

A ht 

A hh 

A ha 

\ 

5n/ R 

A a0 

Axijr 

Axh 

A aa 


5a 

^ J 


= 0 


*80 “ 1 - * a ( x + is e) + rf W 
A e* - -"* 0 * JT + 3 

r K_ 


A 0h = Meh + “2 Q 0h,P 


A 0 a 

ti 

<£ 

p 

*? 

K r 

Q ea,] 

A f 0 


K r 


— M^y 

1 - ? 

2 /“*' 
v H 

A^h 

Kp 

= 3 

K r 


V 

= u vt 

*? 

K r 

as 

= M 0 h 

o 

Qhe,] 





A h\|r 

Kp 

= ~2 V.P 



6\|r 

g0 


Kp 

+ -§ <*i|n|r,P 


k i 

A bh = M hh 


1 - * 2 B (! + *t) 


+ -|(Qhh,P + %h,y) 


A ha = M ha + ^( Q ha,P + Q ha,w) 

kr 


> (B42) 


Ax0 = M 0a + "gT ^aQ.P 
k r 

Aat = ~ m ya. — + 72 Q(X ^ P 

JK r k r 


rs.p 

*ah = M ha + -g(Qah,P + Qah, w) 


A aa ^cax 


- A2 £r( x + ig -) 




+ -|(^io,P + Qaa,w) 


(B4i) 
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In order for this set of equations (eq. (B4l)) to have a nontrivial solu- 
tion, the determinant of the coefficient matrix must vanish. This flutter 
determinant is 


A ee 

A 0l|/ 

A 0h 

A 0a 


A i|n)r 

A i|rh 

H a 

A h0 

A hi|/ 

A hh 

A ha 

A a0 


Aah 

Axa 


= 0 (B43) 


Application of Analysis 

Solution of determinant . - The flutter determinant has been programed for a 
high-speed electronic digital computer. In this program, for a given case, all 
input quantities are fixed except k r , gQ, and A^. A sequence of values of 

reduced frequency k r and propeller structural damping in pitch gQ is used. 
For each combination of values of k r and gQ the flutter determinant is 

expanded to give a fourth degree characteristic polynominal in A^. Points of 
neutral stability of the system are obtained for the combination of values of 
k r and gQ that makes one of the four roots of the characteristic polynominal 

a positive real quantity. The values of ky, gQ, and A^ at this point of 
neutral stability are used herein to compute the flutter frequency ratio 


cd_ = 1 

O) 0 a 


(Bhk) 


the flutter speed ratio 


V _ *>r 1 1 

Ro) 0 R A k r 


(b45) 


and the equivalent viscous damping (see ref. 1) required for neutral stability 


2£q ggA 


(B46) 


Whirl flutter boundaries are known to be quite sensitive to the amount of 
structural damping in the system (see refs. 1, 2, and 5, for example); hence, 
it was not desirable to assume equal damping in all four modes of this analysis. 
This is the reason for the choice of the somewhat lengthy method of solution 
outlined. This procedure allows the wing and propeller-engine damping to be 
varied independently. 
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With the solution of the flutter determinant obtained as outlined here the 
elements in the square matrix in equation (B4l) are completely specified. By 
choosing one of the four generalized coordinates, say q 9 arbitrarily, it is 

possible to solve this set of equations for the ratios of the other generalized 
coordinates to q <x . The relative magnitudes and phases of these complex ampli- 
tude ratios give some indication of the nature of the flutter mode. 

Determination of input quantities .- The input quantities required to solve 
the flutter determinant are listed in equation (b 42) . These quantities were 
obtained as outlined in this section. 

The nondimensional generalized mass parameters are listed in equa- 

tion (B27) . The required integrals were evaluated by using the measured mass 
properties of the model with the computed, uncoupled wing bending and torsion 
mode shapes (fig- 9)- The propeller-engine mode shapes, cp @ and cp^j were 

taken as unity. The wing bending and torsion frequencies employed were those 
calculated for the computed modes. The engine-propeller frequencies, 

and ofy, were uncoupled, measured frequencies. The propeller yaw-to-pitch 
damping ratio, was taken as unity in all cases. The wing structural 

damping coefficient for the first coupled mode was measured and found to be 
0.005* It was assumed that the damping coefficient for both uncoupled wing 
modes was also this value. 

As indicated in table II, there was a small variation in density among the 
data points. In the analysis, a nominal value of p Q = 0.0022 slug/cu ft was 

used. 


The propeller nondimensional generalized aerodynamic forces are listed in 
equation (Bj6) . For no aerodynamic interference between the wing and the pro- 
peller, the following symmetry relations between propeller stability derivatives 
may be assumed: 



Cmq - c n r \ 

S = 



J 


(B47) 


All other propeller derivatives were neglected. The values of the static 
derivatives ^C mQ , ^Zq’ Cm^, C^) used were measured values reported in 
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reference 8 for a windmilling propeller. As in this reference , the values of 
C m ^ employed were those computed by the method of reference 1 6 . These deriva- 
tives are listed in table III for the five propeller blade angle settings used 
in these tests. The values of windmilling-propeller advance ratio, J, and the 
distances from the gimbal axes to the plane of propeller-blade quarter-chord 
points at three-quarter radius, e 0 and e^, are also listed in table III. In 

all cases considered e Q = e^. 

The wing nondimen sional generalized aerodynamic forces are given in equa- 
tion (B39) . The aerodynamic integrals (eq. (b4o)) were evaluated numerically 
by using the calculated bending and torsion mode shapes. The c 7 and a p 

values for the exposed panel were computed from three-dimensional lifting- 
surface theory (kernel function method) for the planform of this model treated 
as a rigid wing at a steady-state angle of attack at zero Mach number. These 
parameters are shown in figure 17- It should be noted that these aerodynamic 
integrals (eq. (b4o)) remain constant throughout the analysis. 
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APPENDIX C 


EFFECT OF GYROSCOPIC COUPLING ON THE FREQUENCIES 
OF VIBRATION MODES 


Derivation of the Frequency Equation 

The effects of gyroscopic coupling on the vibration mode frequencies are 
examined at zero airspeed by determining the coupled vibration frequencies over 
a wide range of propeller rotational speed. The frequency equation for zero 
airspeed is obtained from the flutter determinant (eq. (B43) ) by equating the 
generalized aerodynamic forces Qj_j to zero in the determinant elements 

(eq. (b 42)) and by expanding the resulting determinant. The effect of damping 
on the vibrational frequencies is neglected, such that g is zero in all 
modes. In addition, the factor it/ k r J in the off-diagonal mass terms, which 
is based on the assumption of a windmilling propeller, is put in the form 

•pn 

A. The determinant for the coupled frequencies then becomes: 

b r ^e 


1 - A" 


iM ei(r A 


br^ 


M, 


0h 


M, 


0 a 


-m 


■et 


RO 

b r ^ e 


MJ1 - A 2 i 


<X> C 


• w Rft , 

■iM, A 

ta- bjjCJDQ 


M, 


eh 


Mhhl 1 - A * % 

03 n 


M- 


■ha 


M, 


9 a 


... Rft -v 




ha 


o an 

M 1 - * -f 

“e, 


= o (ci) 


/ \ 0 

where M,- * is defined by equation (B27) . Values of A = — • are desired for 

■ L J CD 

values of SlJ . Expansion of the determinant gives a polynomial in A 2 and 

( — ^ which can be written as: 

\ 9 / 


(SL \ 2 = f M 

\a$ ) \ R / 


2 “Wl 1 ‘ ^ A2 


0 J A i A + a 2 a + a 3 a + A 4 


A" 


+ B 2 A 2 + B, 


(C2) 
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where 


- "^hhHxai 
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^ j V°e ) \“0 j \ cd q 


ai 
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+ *4 + *&?W 


^) 2 + <Ah(^ 


^ M hhHxa *" M ha + ^0h^9aJ%a " ^h^coa, " M 0a^hh 
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1 of hh <xa \u}0 J 1^0 I ^1^0 


“^G^hh^aoa 


avA*" /a^ 
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03 c 


av 


a0^0aVhh(^) ' M fa M hh 


1+1^ 

03q 


B 3 * “f^Hm ' 4^4. + ‘ =^1>WWV + >4“hh " 44* 


Discussion and Application 

Some observations can be made about the behavior of the mechanical system 
from inspection of equation (C2) . For any value of Aj, only one positive value 
of SljcuQ is given by equation (C2) ; thus, the frequency of any one vibration 

mode cannot be the same as the frequency of any other mode at any value of 
ft/o^. Hence, if the frequencies of two vibration modes approach each other as 

varies, the two curves of the frequency plot split and do not cross. 

Such frequency splits are also encountered in other vibrating systems (for 
example, see ref. 17)* The numerator of equation (C2) gives the values of the 
coupled frequencies at 0 /csq of zero. (Note that the yawing mode is uncoupled 

for q of zero.) The denominator of equation (C2) gives three asymptotes: 
A 2 = 0, and the two roots of B^A^ + B 2 A 2 + =0. Thus, with the zeros, 

poles, and uniqueness of for a given value of A, the frequencies of the 

modes can be readily traced for variations of 
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Five cases are considered for illustration: two cases for configuration A, 

and one each for configurations B, C, and D. The results are presented in fig- 
ure 15 for four modes and for only two modes, 9 and i|r. These five cases 
vary widely in the location of the uncoupled engine frequencies relative to the 
wing frequencies. Figure 15 shows that the effects of increasing may 

result in some cases in considerable differences in the values of the coupled 
frequencies from those at H/odq of zero, depending upon the parameters of the 

system. A similar analysis showing the effect of gyroscopic coupling involving 
wing motion for one example is given in reference 18 . 
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TABLE I.- MASS PARAMETERS 


(a) Propeller and power plant (pivoted) 


Configuration 

M p , 

slugs 

*0 or V 

slug-ft2 

Ift, 

slug-ft 2 

s e > 

slug-ft 

A 

o . 2388 

0.0634 

O.OO 858 

0.03801 

B 

.1978 

. 04l48 

.OO 858 

.009128 

C 

•1379 

.01772 

.00858 

-.02849 

D 

•1379 

.01709 

.OO 858 

-.02849 


(b) Nacelle and mounting beam (nonpivoted mass) 


Configuration 


A 

B 

C 

D 


Mass, 

slugs 

Inertia, 

slug-ft 2 

(a) 

Static unbalance, 
slug-ft 
(a) 

0 . l8l4 

0.2237 

-0.1496 

.2224 

.2783 

- . 1962 

.2823 

.3667 

- . 2689 

.2823 

•3674 

- . 2689 


a Values given are relative to the elastic axis of the 


wing. 


TABLE II.- MODEL TEST DATA 





0.0072 


.0074 

I 

.0090 

I 

.0089 


.0091 

.0149 


.0084 


* 

0106 

Jr 

0290 




P> 

slugs/cu ft 


0.0289 


0.0367 


.0292 


Configuration A 


0.0132 


.0091 

>1 

.0108 

l 

.0094 


.0107 

.0112 

.0153 


.0088 


* 

.0120 

J' 

.0324 


0.00230 


* 

.00236 

.00231 

.00230 

.00228 

.00233 


Configuration B 

0.0351 I 0.0299 i 0.00226 


Configuration C 

0.0266 | 0.00233 


Configuration D 


0.0310 


.0199 


0.00233 


rps 


61.5 

41.0 

29.0 

25.0 

21.5 

47.0 

67.5 

48.0 

33.5 

28.0 

23.5 

28.5 

67.5 

50.5 

34.5 

28.0 

24.0 

26.0 

65.5 

42.0 

31.0 
22.8 
68.0 

46.0 

35.0 

31.0 

30.5 

26.0 

63.5 

46.5 

40.0 

36.0 


100 

68.0 

60.5 

53.0 

49.0 


90.5 

72.0 

62.0 

52.0 

50.0 


66.0 

48.5 

37.0 

32.5 

29.0 

27-5 


Vf, 

ft/sec 


125 

121 

126 

135 

148 
138 
138 

143 

144 
151 
162 
154 
143 

149 
149 
151 
168 
l4l 
128 
125 
134 
156 
132 
128 
143 

162 

163 

177 

188 

202 

216 

251 


201 

201 

218 

232 

274 


185 

177 

181 

229 

281 


130 

138 

155 

171 

199 

189 


ff> 

cps 


5-70 

6.58 

7-24 

7.40 

7.67 
6.71 
5.81 
6.66 

7.40 
7-70 
7.89 
8.28 
6.38 

7-30 

7.92 

8.40 

8.46 

8.40 

6.47 
7-78 
8.29 

8.78 
7-06 
8.00 

8.68 

8.79 
9.09 
9-20 
7.60 

8.41 
8.60 
8.94 


6.58 

7-97 

8.46 

8.80 

8.96 


6.16 

7-11 

7-96 

8.44 

8.26 


3-95 

4.69 

5-09 

5.09 

4.86 

5-03 


bo 



TABLE III.- PROPELLER AERODYNAMIC QUANTITIES USED IN THE ANALYSIS 


P 0.75R' 

deg 

J 

e 0 or 

CD 

Cy 

z e 

s 

Cy 

Cm tl 

25 

55 

46 

52 

58 

1-25 

1.80 

2.66 

5-51 

4.21 

-0.51092 

-.31781 

-.32446 

-.32758 

-.33031 

0.0447 

.0382 

.0296 

.0253 

.0213 

-0.232 

-.295 

-.364 

-.402 

-.44-0 

0.0985 

.0964 

.0866 

.0801 

.0734 

0.0996 

.1018 

.1028 

.1032 

.1040 

-0.1574 

-.1126 

-.0738 

-.0525 

-.0361 



















Figure 8.- Sketch of layout of model in tunnel showing measured node lines and frequencies of modes 
that involved primarily wing motion for configuration A , fg - f^ = 9*3 cps. All dimensions in 
inches unless otherwise specified. 
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Figure 11.- Concluded. 



(c) Configuration D. 

Figure 12.- Comparison of theoretical whirl flutter speeds with and without wing flexibility 
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(a) Configuration A, 2£ = 0.029* 

Figure 13 .- Comparison of experimental and theoretical whirl flutter speeds and frequencies. 
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(b) Configuration B, 2t, = 0.032. 
Figure 13 .- Continued. 




(c) Configuration C , 2£ = 0.026. 


Figure 13 •- Continued. 
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Figure l4.- Effect of wing aerodynamics on flutter characteristics for configuration C. 
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Figure l4,- Continued. 
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Figure 15.- Continued. 
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